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1. INTRODUCTION 

Pulsating stars have given important information in classical astronomy and 
are still an active research field, e.g. asteroseismology or helioseismology. By 
studying the rich oscillation stellar spectra, optical astronomy makes predic- 
tions for the stellar equations of state (EOS), the evolutionary age of the stars, 
and in the case of Cepheids one can make exact predictions for the distance of 
the object. 

For relativistic astrophysics the pulsations of relativistic stars are of great 
importance since through them one can discuss stability properties of compact 
objects and, moreover, pulsating stars can be a promising source of gravita- 
tional waves. It is expected that during the first few seconds which will follow a 
supernova explosion the newly formed neutron star will pulsate wildly and this 
pulsation will be mainly damped due to the emission of gravitational radiation 
which will carry away the characteristic signature of the collapsed object. The 
energy stored in the pulsation will be of the same order as the kinetic energy 
of the collapse and in this way a significant part of the original mass-energy of 
the newly formed neutron star will be radiated away as gravitational radiation. 
With the present sensitivity of the resonant detectors it is quite sure that the 
gravitational waves will be detected if a supernovae explosion takes place in 



(*) Review talk for the proceedings of the Les Houches School on Astrophysical 
Sources of Gravitational Waves Eds J-A. Marck and J. P. Lasota, Springer- Verlag 
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the Local Group of galaxies. Depending on the details of the collapse EJ, with 
the future generation of resonant detectors it is possible that events up to the 
distance of the Virgo cluster to be detectable. 

The study of stellar pulsations in General Relativity (GR) has been initiated 
nearly thirty years ago by Thorne and his collaborators [g, IJ [|, [|, [f|, and 
was intended mainly in the calculation of the / (fundamental) mode since this 
is the mode through which most of the gravitational radiation of the star is 
radiated away. All the other fluid modes, g (gravity), p (pressure), s (shear), 
t (toroidal), i (interface) modes, can be calculated with quite high accuracy 
with the Newtonian dynamics since they don't emit significant amounts of 
gravitational radiation. This Newtonian like picture has changed significantly 
in the last ten years due to some new ideas by Kokkotas and Schutz 0, || and 
Chandarsekhar and Ferrari || but also due to subsequent work by Andersson 
and Kojima. The new approach has been based on the fact that the stellar 
pulsations in GR should not be treated in a semi-Newtonian manner because 
then some important features of the theory will be overlooked since they are 
related uniquely to the relativistic nature of the problem and they do not exist 
in the Newtonian theory. These unique new features are: 

• the existence of a new set of modes, the gravitational wave modes (w- 
modes), which are purely spacetime modes and which cannot be seen in 
Newtonian theory g | 0, [H], @ [l| || ||. 

• that the stellar pulsations can be treated as a scattering problem [^) i.e. in 
the perturbation equations only the variations of the gravitational field 
will enter while the fluid perturbations will be absent. This approach 
exists even in the limit of Newtonian stars (l6[ [l7| . 

• that in the slow rotation limit the polar modes excite the axial ones due 



to the dragging of the inertial frames |18, [L9| [20| [21) 

• that the axial modes which were thought not to exist for non-rotating 
stars do exist ^2) and that they show a spectrum similar to the w-modes 
|2j| and that, more or less, they can be considered on an equal footing 
as the polar modes |l3). Recently, it has been shown [^4| that these 
modes can be excited and emit as much energy as an oscillating black- 
hole. 

These new developments in the last decade in the theory of relativistic stellar 
pulsations will be the subject of the following sections. We will not discuss 
stability since for non-rotating stars there is no significant progress in the last 
ten years and thus the reader can refer to a review article by Schutz [p5[. 



2. FLUID MODES 



The fluid modes are modes that exist in both Newtonian and GR treatment 
of the stellar pulsation and their main characteristics do not depend on the 
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specific gravitational theory. If we consider an unperturbed fluid it will be 
described in both theories by the mass and momentum conservation equations 
and an equation which describes the gravitational field (Poisson or Einstein 
equations). By assuming a small variation from sphericity we can practically 
describe all fluid and dynamical quantities in terms of spherical harmonics while 
a harmonic dependence on time can also be used. In Newtonian theory all the 
quantities which describe the variations of both the fluid and the potential are 
scalar quantities and their decomposition in spherical harmonics has only one 
parity. In GR, the variations of the spacetime metric are tensor quantities and 
thus their decomposition in spherical harmonics has two parities, odd and even 
ones. The even parity or polar perturbations are similar to the Newtonian ones 
and describe spherical deformations while the odd parity or axial ones (called 
also toroidal) are degenerate in Newtonian theory. They can set the star into 
a continuous non-varying differential rotation. Nevertheless, this degeneracy is 
erased in the presence of rotation, magnetic fields, or non-zero shear modulus. 
Thus, in principle, in the absence of rotation in both theories these modes are 
unimportant because they cannot describe time dependent motions and as a 
result there is no emission of gravitational radiation. We shall discuss these 
modes in the next section where we shall see that the Newtonian approach to 
the problem is not complete. The perturbation equations have the general form 



where W(r) is a function of the radial, £ r , and tangential, £g, components 
of the fluid displacement vector. Z(r) is a function of the components which 
describe the variation of the gravitational field [(^ £ is the harmonic index (m 
is not present because we have a (2£ + l)-fold degeneracy due to the rotational 
symmetry of the equilibrium structure around any arbitrary axis) and u) is 
the pulsational frequency. 2?i and T>i are both wave operators in GR with 
propagation velocities, the acoustic one for the first and the speed of light for 
the second equation (in Newtonian theory T>2 is not a wave operator thus we 
have no emission of gravitational waves). In GR outside the star the equation 
ph reduces to either Zerilli (polar) or Regge- Wheeler (axial) wave equation 
|2S| ] depending on the parity. The equations (jl]) and (|2|) together with a set 
of appropriate boundary conditions formulate an eigenvalue problem with u> 
being the eigenvalue and W(r) and Z{r) the corresponding eigenfunctions. 
The boundary condition for the existence of quasi-normal modes (QNM) is that 

( ) In the Newtonian theory they are the variation of the potential 5U and its 
derivative, while in GR they are components of the perturbed part of the metric, e.g. 
for a certain choice of the gauge, htt and hgg . In both cases the system of differential 
equations is of 4th order pq, in the gauge used in |H| the system is of 5th order but 
it has been proven |?7j that by a gauge transformation it can always reduced to a 4th 
order one. A gauge invariant presentation by Moncrief 0] has never been used up 



2W + F (W, Z', W, Z- £,uj) = 



(1) 



V 2 Z + G(W',Z',W,Z;£,lu) = 



(2) 



to now. 
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at infinity there is no incoming gravitation radiation. In this way a discrete 
spectrum of frequencies for a certain stellar model can be revealed. 

The richness of the spectrum depends on the complexity of the stellar model 
and each new feature that we add to the stellar model results in a new family 
of modes, [p] The simplest possible model for a neutron star is a non-rotating 
fluid ball at zero temperature. If the fluid has constant density this model 
supports only one pulsation mode for each multipole i. This mode, which in 
Newtonian theory has oscillation frequency 

lo(R 3 /M) 1/2 = y/2£(£~~l)/(2£+l) w 0.894 for I = 2 , (3) 

where R and M are the radius and mass of the star in units where c = G = 



1, was first studied by Kelvin 31 and is usually referred to as the /-mode 
(fundamental mode). It is distinguishable by the fact that there are no nodes 
in the corresponding cigenfunctions inside the star. In a way, the /-mode is 
due to the interface between the star and its surroundings. The eigenfunctions 
of such modes would typically have maxima at the interface and fall off away 
from it. This is exactly the character of the /-mode: It reaches maximum 
amplitude at the surface of the star (see for example Figures 7-8 in ||). A 
comprehensive set of quadrupole /-mode calculations has been carried out for 
13 equations of state by Lindblom and Detweiler (3^| . The frequencies are in 
the range of 1.5-3.5 kHz and the damping times are of the order of 0.1-0.5 sees, 
i.e. the star undergoes a few hundred oscillations before it damps out. 

A somewhat more realistic star consists of a perfect fluid. Then one must 
specify the equation of state, and most studies to date have been for simple 
polytropes. For this stellar model a second set of modes - the p-modes, the 
restoring force of which is pressure - exists. The oscillation frequencies of the 
p-modes depend directly on the travel time for an acoustic wave across the 
star. This would typically lead to higher oscillation frequencies starting from 
5-6 kHz while their damping times are significantly longer than those of the 
/-modes and increase with the order of the mode. 

When the temperature of the star is non-zero a further set of modes comes 
into play. These modes are mainly restored by gravity and are consequently 
referred to as g- modes. For a star in convection, i.e., when the entropy is 
constant, the (/-modes are all degenerate at zero frequency. In general, however, 
their oscillation frequencies depend directly on the central temperature of the 
star and are typically smaller than a few hundred Hz. 

The three families of fluid modes discussed so far, the /-, p- and g-modes, 
all belong to the class of polar modes. For these models there are no fluid axial 
modes, since the stellar models discussed are all somewhat unrealistic. If even 
more realistic stellar models are considered, for example it is expected that 
neutron stars will have a kilometer-thick crust which will crystallize, then, if 
the shear modulus in the crust is non-zero, axial modes do exist |33L p3. There 



( 2 ) The naming of the various families of fluid modes comes from a monumental 
work by Cowling in early '40s fjcjj . 
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will be also families of modes directly associated with the various interfaces 
inside the star which will not be discussed here and for a detailed discussion 
the reader should refer to an excellent presentation in j3f| . 

The maximum pulsation energy is stored in the /-mode on which the fluid 
parameters undergo the largest changes, typical values are of the order of 10 53 — 
10 54 ergs |32|]. If the gravitational radiation is the only damping mechanism 
of the pulsations (which is approximately true ]35[|) then for the Newtonian 
ones the emission rate can be calculated by the quadrupole formula ||_|36| 
and the damping rate will be just E/E (for the analytic formula sec ]37|). 
In GR approach the pulsation frequency is complex and the imaginary part 
corresponds to the damping of the pulsation. The calculations of the /-mode 
frequencies for the same mean density stars have shown that both Newtonian 
and GR predict more or less the same values. This picture is also true for all 
the other families of fluid modes and thus they can be studied in Newtonian 
theory |35|]. The previous picture is partially true for the damping rate since 
when the star is fully relativistic {e.g. M 1.4M© and R m lOKm) the 
Newtonian theory predicts 2-3 times faster damping of the oscillations |38| . The 
post-Newtonian (pN) treatment of the problem reduces this gap significantly 
|3S| , fh| |d| and it is a promising approach especially for the study of rotating 
stars. 

As it has been seen by equation (^|) the /-mode frequency scales with the 
mean density of the star, while its damping time (due to radiation reaction) 
depends on the compactness of the star, i.e. the more compact the star is the 
faster the oscillations damp out (for a discussion of the damping times including 
viscosity see Q ) . These two properties of the /-mode frequency are useful for 
the estimation of the masses and radii of pulsating neutron stars. 

We shall conclude this section by mentioning a very useful approximation 
for the study of the fluid modes, the so called Cowling Approximation in which 
the perturbations of the potential (or the metric) are set to zero and thus the 
whole problem is described with just one wave equation for the fluid. In this 
approximation both frequencies and damping times can be revealed with an 
error usually less than 10% |l3], |44|, An additional advantage is that the 
simple wave equation can be studied even analytically in the WKB approxi- 
mation and the properties of the various families of modes described above can 



be visualized in an elegant way 46 



3. SPACETIME OR W-MODES 



As it has been shown in the previous section the treatments of stellar pulsations 
in Newtonian theory and in GR do not differ significantly. Actually, when the 
pN correction were considered even the gap in the estimation of the damping 
times has been considerably reduced. Thus the natural question is what is 
actually the role of the equation (||) being a wave one in GR? Is it just for the 
direct calculation of the damping times during the solution of the eigenvalue 
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problem and not afterwards from the ratio E/E? The answer is clearly no. 
Spacetime has it own dynamics and does not play a passive role as just the 
medium in which gravitational waves propagate; moreover, it has its own unique 
spectrum. 

Let us start by assuming a simplified version of the problem. We shall con- 
sider the Inverse of the Cowling Approximation (ICA) JlJ] which has been 
discussed earlier, i.e. instead of freezing the spacetime perturbations we shall 
freeze the fluid ones. Then the whole problem is described by just one second 
order wave equation both inside and outside the star (as for the axial pertur- 
bations), which together with the appropriate boundary conditions (regularity 
of the perturbed quantities in the center and only outgoing waves at infinity) 
forms a well posed eigenvalue problem. The spectrum should have, in princi- 
ple, a lot of similarities with that of QNM of black holes [^Tj, but since the 
boundary conditions at the center of the star are different than those at the 
horizon for a black-hole, the two spectra should not be identical. Nevertheless, 
the characteristic times of the perturbation described by the wave equation in 
the ICA will be proportional to the wave speed which in this case is the speed 
of gravitational waves and thus the perturbations should be of higher frequency 
and they should be damped out very fast (we remind you that for the fluid pul- 
sation the characteristic time was related to the speed of the acoustic waves). 
In this way we reveal a new family of modes which is related directly to the 
spacetime and which can never be seen through the Newtonian approach to 
the problem. 

It took us nearly 20 years (from the original work by Thorne) to recognize 
this natural GR extension of the stellar pulsation spectrum. This approach 
has initiated ten years ago when Kokkotas and Schutz Q| have studied the 
QNMs of a simple "toy" model consisting of two strings, one finite and one 
semi-infinite coupled together with a spring; this system mimics very well the 
star-spacetime system. They have found that there exist two families of modes, 
one with slow damping (like the fluid modes) which were mainly modes of the 
finite string (star) and another family with very fast damping which was mainly 
modes of the semi-infinite string (spacetime). The damping time of the QNMs 
of this second family was shorter as the coupling (the spacetime curvature) 
of the two strings were looser. Such instructive "toy" models can now easily 
be constructed and all of them suggest the existence of this second family of 
modes |f,[||f(J. 

A few years later this new family of spacetime modes has been found by the 
same authors Q . They named these new modes "gravitational wave" modes or 
w-modes. [PJ The characteristic properties of these modes, as it was expected 
by the previous analysis and the "toy" model, were: (a) high frequencies (8-12 
kHz), (b) fast damping times (0.02-0.1 msecs) which decrease with the order 
of the mode and/or with the compactness of the star (as the "toy" model 

( 3 ) It should be mentioned that Kojima ju^] had earlier found a few of the w-modes 
but he had not revealed the whole spectrum and its properties. 
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Fig. 1. — A graph which shows all the ui-modes: curvature, trapped and interface 
both for axial and polar perturbations for a very compact uniform density star with 
M/R = 0.44. The black-hole spectrum is also drawn for comparison. As the star be- 
comes less compact the number of trapped modes decreases and for a typical neutron 
star (M/R — 0.2) they disappear. The Im(oj) = 1/damping of the curvature modes 
increases with decreasing compactness and for a typical neutron star the first curva- 
ture mode nearly coincides with the fundamental black-hole mode. The behaviour of 
the interface modes slightly changes with varying compactness (cf. Figure 3). The 
similarity of the axial and polar spectra is apparent. 



predicted), and (c) practically no significant fluid motions inside the star. 

The w-modes exist for both polar and axial perturbations since they do not 
depend on the perturbations of the fluid, actually in the black-hole perturba- 
tions the spectra in the two cases coincide. Nevertheless, for stars the boundary 
conditions near the center are different than those for black-holes (as well as 
the form of the potential inside the star) and naturally we should not expect 
coincidence, but instead a similarity of the axial and polar u>-mode spectra as 
it has been shown in |f4j , see also Figure [l]. So in what follows we shall 

not discriminate the axial and polar perturbations. 

A way to understand the nature of these modes is the "trapping" of the 
gravitational waves by the spacetime curvature inside the star, ft is easy to 
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see how this may happen if one plots the gravitational- wave speed as measured 
by an observer at infinity: g u as a function of r. This has a minimum at the 
center of the star, and the interior w-modes would be trapped in this "bowl" 
of curvature. Moreover, such modes would naturally be concentrated at the 
center of the star, which agrees well with the eigenf unctions constructed in |lq] 
(see also an extensive discussion in |l3|). Hence, it makes sense to refer to these 
modes as " curvature modes". 

A naive but useful argument leads to standing wave solutions inside the star, 
essentially of the form sin(wr). If this were the true form of the solutions, and 
the modes only leak out slowly through the surface, we would get 

Lo n R = niT , n = l,2, ... (4) 

where R is the stellar radius. This argument is far too simplistic, but it is 
rewarding to find that two of its predictions agree well with the results for the 
curvature modes: There would exist an infinite sequence of such modes, and 
their pulsation frequencies Re(w„i?) would be separated by 7r (cf. Figure 5 in 



13 ). This dependence on the size of the star was also predicted by the simple 
"toy" model in 0. 

When the star is made increasingly compact, a few of the curvature modes 
clearly change character and should be considered as trapped in the potential 
well that arises inside the black- hole potential barrier (R < 3M). In this 
regime even the fluid /-mode should also be considered as a trapped mode 
and shares the properties of the w-modes (cf. Figure 1 in [|l3)). It will be 
instructive to name all the modes which are trapped because of the potential 
well trapped spacetime modes. Their number is finite and depends on the 
depth of the potential well. The existence of trapped modes has been first 
discussed by Detweiler |[l| but the interest on them has been revived with the 
work of Chandrasekhar and Ferrari |^2| for the axial stellar perturbations. The 
existence of trapped modes for both polar and axial perturbations has been only 
recently shown jl2[ |l3|, [l4| but they are unlikely to be of any great astrophysical 
relevance, such compact stars will probably never form. 

The other family of w-modes discovered by Leins et al. |ll| was not pre- 
dicted before by any "toy" model, and it is characterized by the extremely fast 
damping and the quite low frequencies; they named them u>//-modes. The ex- 
istence of this family of w-modes may be more directly due to the discontinuity 
at the surface of the star. Then the mode-eigenfunctions need not be local- 
ized in the star. Rather, these modes would be similar to modes for acoustic 
waves scattered of a hard sphere. One would typically expect such modes to be 
short-lived compared to modes trapped inside the star. Only a finite number of 
modes exist for each multipole £, and there may be purely imaginary ones. The 
latter feature suggests that one should perhaps not be surprised to find stellar 
modes "emerge" from points on the imaginary-frequency axis as in Figure [| 
here. This evidence seems to be compelling and makes the association of the 
second family of w-modes and the interface at the stellar surface likely Hence, 
it will be instructive to call these modes interface modes. In principle, the 
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Fig. 2. — (left) Re u)(R z /M) 1/2 and (right) Im lj{R 3 /M) 1/2 as functions of the 
compactness of the star for the /-mode (solid and denoted by f), and the first few 
ui-modes for polar (solid) and axial (dashed) perturbations. Notice the beautiful 
example of avoided crossings in the pulsation frequencies. 



conjecture that only a finite number of these modes exist should be testable. 
But at present numerical difficulties restrict calculations to lm(u>M) < 1.25 or 
so. Much better numerical schemes, or other formulations of the problem, are 
required to test this prediction. 

Finally, we shall discuss some of the new results induced by the slow rota- 
tion of neutron stars. As we have mentioned earlier the eigenfrequencies are 
degenerate with respect to m at fixed t. In the presence of rotation, however, 
the degeneracy with respect to m is removed and different modes are mixed 
with each other. The frequency and damping times of the co-rotating /-mode 
increase with the stellar angular velocity while those of the counter-rotating 
mode decrease |2C|. Moreover, the axial mode with £,m is coupled with the 
polar mode with £ ± 1, m and vice versa jl8], [l9| |2lJ. This results from the 
dragging of the inertial frames (Lense-Thirring effect) by the rotation of the 
star. 



4. EXCITATION AND PARAMETER ESTIMATION 

In the previous sections we have seen all the new developments in our under- 
standing of the role of the spacetime for pulsating stars. Nevertheless, many 
important questions remain, and the most important of all of these is whether 
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Real part 



Fig. 3. — \m(uoM) vs Ke(uM) for the polar (solid) and axial (dashed) interface 
modes. Only the first of these modes exist for all values of stellar compactness. 
Further modes arise for sufficiently compact stars. The numerical calculation becomes 
difficult for large Im(tjM) and consequently we have only partial data for the third 
of the polar modes modes. Arrows indicate the direction of increasing compactness. 



the w-modes can contribute to observable gravitational waves, and thus play 
a role in astrophysics. To test this, Andersson and Kokkotas |24j have studied 
scattering of axial gravitational wave-packets by a compact star. The problem 
is similar to the black hole one that was studied by Vishveshwara in early '70s 



52 . The axial problem corresponds to a single wave equation with an effective 
potential. It is remarkable that what must be considered a basic exercise in 
numerical analysis can provide us with new and important information. 

The result of a typical simulation is shown in Figure ^. The exponential 
ring-down at late times (from t k 100M) corresponds to the first axial w-mode 
(Figure ^B-D). A theoretical waveform based on the slowest damped axial w- 
mode for this star fits the numerical one perfectly. After the QNM ring-down 
(at very late times) the waveform is dominated by a power-law tail. This is 
exactly what one would expect ]53[ | . It is important to be mentioned that the 
last part of the signal would be the same if the star was replaced by a black 
hole of the same mass. 

Similar simulations for ultracompact stars (R < 3M) have been also per- 
formed (Figure ||A-C), for these there will be a few trapped modes. In this case 
we generally find that the first mode "above the peak" [ with (ujM) 2 > V max , 
this is not a trapped mode] dominates the radiation initially. Most of the en- 
ergy goes into this mode, but since it is more rapidly damped than the trapped 
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Fig. 4. — The response of two uniform density stars (Al/R = 0.44 left and Al/R — 
0.20 right) to a Gaussian pulse of gravitational waves. The top shows the actual 
form of the axial perturbation as seen by a distant observer, while the lower panel 
shows the same function on a logarithmic scale. The QNM ringing that is apparent 
after t « 100M (right) corresponds to the first axial curvature w-mode, while there 
are no trapped lo-modes. The power-law tail dominates after t « 150M. In the left 
picture we can see a series of QNMs which correspond both to curvature and trapped 
w-modes. 



ones, the latter dominate the radiation at later times. But as we mentioned 
earlier the ultracompact stars are interesting only from a theoretical point of 
view: We can learn little astrophysics from them. An extension of the above 
work is to study the excitation of both fluid and spacetime modes, i.e. polar 
perturbations, there one can exactly determine the amount of energy radiated 
due to fluid (/, p-modes) and spacetime modes ]5r| . 

( 4 ) Recently, Borelli and Ferrari p3| came to similar results about the excitation 
of the trapped axial modes by a mass falling onto a compact star. Their method is 
similar to the one employed by Kojima |^5| for the study of stellar resonances due to 
orbiting particles. 
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The simulations suggest that the w-modes will be relevant in many dynamical 
processes involving neutron stars. The modes ought to be excited when a 
stellar core collapses to form a neutron star. Much of the initial deformation 
of spacetime could then be radiated away in terms of u>-modes. The modes 
should also be excited when two neutron stars collide or at the final stages of 
binary neutron star coalescence. 

It seems clear that the w-modes can play a role in many astrophysical sce- 
narios, but will we be able to observe them with future gravitational- wave 
detectors? This is a key question, the answer of which demands much more 
detailed calculations than the ones presented here and first of all 3D fully GR 
codes for gravitational collapse or fully GR evolution codes for the binary neu- 
tron star coalescence. But let us nevertheless suggest a handwaving argument: 
A spectral analysis (see Figure |^) shows that the w-modes will be excited to 
roughly the same level as the modes of a black hole would be in a similar situa- 
tion. The black-hole modes are expected to be detectable pl\ , so the situation 
looks promising also in the case of stars. 




5 10 15 20 25 30 35 

frequency ( kHz ) 



Fig. 5. — The power spectrum for the simulation in Figure [l|B-D. The axial power 
spectrum (solid line) is compared to the corresponding one for a black hole (dashed 
line). We also show the power in the initial Gaussian pulse (another solid line) that 
has been used to normalize the other spectra. Note that the star mode is excited 
to roughly the same level as the black-hole one. A significant part of the initial 
energy clearly goes into quasinormal-mode ringing. The insertion shows the axial 
power spectrum in logarithmic scale. From this it is clear that the first three w- 
modes are excited (the corresponding frequencies are uiM « 0.35, 0.82 and 1.25 
which corresponds to ui w 8.36, 19.6 and 29.8 kHz) 
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In fact, the frequencies of the /-modes (around 1-2 kHz) make them well 
suited for detection by the operating resonant bar detectors and the recently 
proposed spherical solid-mass detectors, e.g. TIGA |33], The detection of 
the w-modes (and even of p-modes) which have frequencies of the order of 8-12 
kHz will be possible by the proposed arrays of small resonant detectors p5| . 
If a substantial fraction of the binding energy of a neutron star were released 
in ii;-modes then these detectors could well see such events out to the Virgo 
cluster. 

Suppose that we detect a gravitational-wave signal from a compact star, 
what can we hope to learn from it? In what follows we shall suggest a possible 
scheme that looks very promising. The idea is the following: Assume that we 
detect a signal and manage to extract both the fundamental polar w-mode and 
the fluid /-mode from it. The discussion in the previous two sections suggests 
that: 

• the oscillation frequency of the /-mode scales with the mean density of 
the star ^M/R 3 (cf. Figure g-I), 

• the damping rate of the w-mode scales with the compactness ratio of the 
star M/R [| (cf. Figure |-IV). 

These properties are illustrated in Figure ^. In principle, one should therefore 
be able to infer both the mass and the radius of the star from the observed 
data. Additionally, the information from the damping time of the /-mode 
(Figure §TI) and the frequency of a w-mode (Figure ||-III) will pose strict 
constraints on the possible neutron star EOS. We are not aware of any other 
scheme that would enable us to extract such a precise information, and put 
as strong constraints on the nuclear equation of state. Moreover, if the star 
is rotating the spectrum changes and the degree of deviation from the non- 
rotational case could help us to estimate the rotation rate of the star just as 
for black- holes |p7| . 

The idea seems simple enough, but will it be useful in practice? We have 
constructed a number of different neutron star models with realistic EOS from 



those listed in 1 32 and we have determined the /-mode and the slowest damped 
polar w-mode for each of these models. The relevant data are graphed in Figure 
|[ The theoretical spectra are well fitted by the following two relations: The 
oscillation frequency varies with the average density of the star as 



uj(kHz) = 0.39 + 44.45 ( ) (5) 



while the damping rate of the w-mode behaves as 



M fM\ 2 



T{msec) = 0.1 - — + 2.96 — (6) 
R V R J 

(M, R in kms). If the two relations (|j) and (^) provide reasonably accurate 
estimates for both M and R for all stars in our dataset, then this idea passes 
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Fig. 6. — (I) The pulsation frequency of the /-mode as a function of the mean 
density of the star, (II) the damping of the same /-mode as function of the stellar 
compactness, (III) the frequency of the slowest damped w-mode as function of the 
stellar compactness and (IV) the damping rate of the same w-mode as a function of 
the stellar compactness. The symbols stand for the corresponding EOS i.e. A for 
Pandharipandi (1971), C for Bethe and Johnson (1974), E for Moszkowski (1974), F 
for Arponen (1972) and G for Canute and Chitre (1974). 



a first test. The procedure passes this simple test with flying colours: The 
mass and the radius of each star can be determined to within 5% by inverting 
(|) and (§). The results of such a test are shown in Table |, where we have 
tested whether the above relations can estimate the masses and radii of various 
polytropes. 

The evidence provided here is, of course, only an indication that this idea 
can be useful in a real detection situation. To investigate this in more detail 
one must incorporate the estimated effect of statistical errors and measurement 
ones. It will, for example, be much more difficult to infer the w-mode damping 
rate from a data set than to find the /-mode pulsation frequency (the observ- 
ability of a periodic signal buried in noise scales roughly as the square-root of 
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Table I. — Test of the estimation of parameters hypothesis on polytropic models 
with N=0.8, 1 and 1.2, the values in parentheses are the estimated values from the 
relations (||) and (^j). 



N 


R (Km) 


M/Mq 


M/R 


L0 f (kHz) 


t w (msecs) 


0.8 


10.21 (10.54) 


0.92 (0.97) 


0.133 (0.135) 


1.933 


0.0152 


0.8 


9.49 (9.61) 


1.36 (1.35) 


0.211 (0.207) 


2.532 


0.0238 


1 


8.86 (8.35) 


1.27 (1.16) 


0.206 (0.211) 


2.870 


0.0236 


1 


7.42 (7.28) 


1.35 (1.31) 


0.266 (0.269) 


3.668 


0.0493 


1.2 


10.48 (9.89) 


1.44 (1.48) 


0.203 (0.221) 


2.540 


0.0276 


1.2 


8.97 (9.45) 


1.46 (1.45) 


0.240 (0.254) 


3.128 


0.0422 



the number of observed cycles), and also to obtain fits similar to (g|) and (g) for 
all known realistic EOS. If such relations prove as robust as the ones we have 
obtained here then the suggested scheme looks truly promising. In principle it 
will be possible even to identify exotic type of stars like the boson stars which 
also show a w-mode spectrum |56) . 

5. CONCLUSIONS 

In this review we have shown that the study of the stellar pulsations as sources 
of gravitational waves should be performed in GR since in the Newtonian grav- 
ity apart from deviations in the strong field regime, a lot of useful physical 
information is overlooked. 

The existence of new families of modes (not seen in Newtonian theory) which 
are not just another branch of the possible mathematical solutions of the eigen- 
value problem but instead are modes which can be excited and detected enrich 
the information that we can get from the signals of pulsating stars. 

The assertions presented here must be tested by more detailed, fully general 
relativistic simulations. This is, in fact, an interesting challenge for numerical 
relativity. A relativistic description of gravitational collapse to form a neutron 
star, or the merger of two stars, should tell us whether the iu-modes are of 
observational relevance or not. Given the present importance for TIGA, fully 
relativistic simulations are urgently required. 

Moreover, there is still a lot of work to be done on the understanding of 
stellar pulsations. For example, we have not discussed the rotating case, which 
in GR is tractable only in the limit of slow rotation. Thus a natural question is: 
do the w-modes induce new instabilities or enchance the classical gravitational 
radiation instability for fast rotating stars? The answer to this question is of 
great importance for the gravitational wave astronomy and, due to the lack of 
full GR rotating star solutions, one should try to attack the problem via pN 
approach |67|| . 
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